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Based on a calculation of neural decoherence rates, we ar- 
gue that that the degrees of freedom of the human brain that 
relate to cognitive processes should be thought of as a classical 
rather than quantum system, i.e., that there is nothing funda- 
mentally wrong with the current classical approach to neural 
network simulations. We find that the decoherence timescales 
(~ 1CT 13 - 10" 20 seconds) are typically much shorter than 
the relevant dynamical timescales (~ 10 -3 — 10 _1 seconds), 
both for regular neuron firing and for kink-like polarization 
excitations in microtubules. This conclusion disagrees with 
suggestions by Penrose and others that the brain acts as a 
quantum computer, and that quantum coherence is related 
to consciousness in a fundamental way. 



I. INTRODUCTION 

In most current mainstream biophysics research on 
cognitive processes, the brain is modeled as a neural net- 
work obeying classical physics. In contrast, Penrose 
and others have argued that quantum mechanics may 
play an essential role, and that successful brain simula- 
tions can only be performed with a quantum computer. 
The main purpose of this paper is to address this issue 
with quantitative decoherence calculations. 

The field of artificial neural networks (for an introduc- 
tion, see, e.g., is currently booming, driven by a 
broad range of applications and improved computing re- 
sources. Although the popular neurological models come 
in various levels of abstraction, none involve effects of 
quantum coherence in any fundamental way. Encouraged 
by successes in modeling memory, learning, visual pro- 
cessing, etc. many workers in the field have boldly 
conjectured that a sufficiently complex neural network 
could in principle perform all cognitive processes that we 
associate with consciousness. 

On the other hand, many authors have argued that 
consciousness can only be understood as a quantum ef- 
fect. For instance, Wigner Q suggested that conscious- 
ness was linked to the quantum measurement problemFL 
and this idea has been greatly elaborated by Stapp M. 
There have been numerous suggestions that conscious- 
ness is a macroquantum effect, involving superconduc- 



tivity fll^l , superfluidi ty Jl3| ], electromagnetic fields p4[ , 
Bose condensation |]l5|,|l6| , superflourescence [ fl7f or some 
other mechanism [ |18|J_19| . Perhaps the most concrete one 
is that of Penrose proposing that this takes place 
in microtubules, the ubiquitous hollow cylinders that 
among other things help cells maintain their shapes. It 
has been argued that microtubules can process informa- 
tion like a cellular automaton [^0) , and Penrose suggests 
that they operate as a quantum computer. This idea has 
been further elaborated employing string theory methods 

The make-or-break issue for all these quantum mod- 
els is whether the relevant degrees of freedom of the 
brain can be sufficiently isolated to retain their quan- 
tum coherence, and opinions are divided. For instance, 
Stapp has argued that interaction with the environment 
is probably small enough to be unimportant for cer- 
tain neural processe s |28[| , whereas Zeh j2j|, Zurek pfj| , 
Scott plj , Hawking [p2| and Hepp |Q have conjectured 
that environment-induced coherence will rapidly destroy 
macrosuperpositions in the brain. It is therefore timely 
to try to settle the issue with detailed calculations of the 
relevant decoherence rates. This is the purpose of the 
present work. 

The rest of this paper is organized as follows. In Sec- 
tion ||, we briefly review the open system quantum me- 
chanics necessary for our calculations, and introduce a 
decomposition into three subsystems to place the prob- 
lem in its proper context. In Section III, we evaluate 
decoherence rates both for neuron firing and for the mi- 
crotubule processes proposed by Penrose et al., relegating 
some technical details to the Appendix. We conclude in 
Section IV by discussing the implications of our results, 
both for modeling cognitive brain processes and for in- 
corporating them into a quantum-mechanical treatment 
of the rest of the world. 



II. SYSTEMS AND SUBSYSTEMS 

In this section, we review those aspects of quantum 
mechanics for open systems that are needed for our cal- 
culations, and introduce a classification scheme and a 
subsystem decomposition to place the problem at hand 
in its appropriate context. 



1 Interestingly, Wigner changed his mind and gave up this 
idea jl^] after he became aware in of the first paper on deco- 
herence in 1970 lllf. 
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A. Notation 

Let us first briefly review the quantum mechanics of 
subsystems. The state of an arbitrary quantum system 
is described by its density matrix p, which left in isolation 
will evolve in time according to the Schrodinger equation 



(1) 



It is often useful to view a system as composed of two 
subsystems, so that some of the degrees of freedom cor- 
respond to the I s * and the rest to the 2 nd . The state of 
subsystem i is described by the reduced density matrix 
Pi obtained by tracing (marginalizing) over the degrees 
of freedom of the other: p\ = ti2P, P2 = trip. Let us 
decompose the Hamiltonian as 



H = Hi + H2 



Hi, 



(2) 



where the operator Hi affects only the 1 st subsystem 
and H2 affects only the 2 nd subsystem. The interaction 
Hamiltonian Hint is the remaining nonseparable part, de- 
fined as Hi n t = H — Hi — H2, so such a decomposition 
is always possible, although it is generally only useful if 
Hi n t is in some sense small. 

If H lnt = 0, i.e., if there is no interaction between 
the two subsystems, then it is easy to show that pi — 
—i[Hi, pi]/h, i — 1,2, that is, we can treat each subsys- 
tem as if the rest of the Universe did not exist, ignoring 
any correlations with the other subsystem that may have 
been present in the full non-separable density matrix p. 
It is of course this property that makes density matrices 
so useful in the first place, and that led von Neumann 
to invent them [|4|: the full system is assumed to obey 
equation (]]]) simply because its interactions with the rest 
of the Universe are negligible. 



B. Fluctuation, dissipation, communication and 
decoherence 

In practice, the interaction H lnt between subsystems 
is usually not zero. This has a number of qualitatively 
different effects: 

1. Fluctuation 

2. Dissipation 

3. Communication 

4. Decoherence 

The first two involve transfer of energy between the sub- 
systems, whereas the last two involve exchange of infor- 
mation. The first three occur in classical physics as well 
- only the last one is a purely quantum-mechanical phe- 
nomenon. 

For example, consider a tiny colloid grain (subsystem 
1) in a jar of water (subsystem 2). Collisions with water 



molecules will cause fluctuations in the center-of-mass 
position of the colloid (brownian motion). If its initial ve- 
locity is high, dissipation (friction) will slow it down to 
a mean speed corresponding to thermal equilibrium with 
the water. The dissipation timescale Tdi ss , defined as the 
time it would take to lose half of the initial excess energy, 
will in this case be of order r co n x (M/m), where r co n is 
the mean-free time between collisions, M the colloid mass 
M and m is the mass of a water molecule. We will define 
communication as exchange of information. The infor- 
mation that the two subsystems have about each other, 
measured in bits, is 



Si 



s, 



(3) 



where Si = —tiiPi log pi is the entropy of the i th subsys- 
tem, S = — trplogp is the entropy of the total system, 
and the logarithms are base 2. If this mutual informa- 
tion is zero, then the states of the two systems are un- 
corrected and independent, with the density matrix of 
the separable form p = pi ® P2- If the subsystems start 
out independent, any interaction will at least initially 
increase the subsystem entropies Si, thereby increasing 
the mutual information, since the entropy S of the total 
system always remains constant. 

This apparent entropy increase of subsystems, which 
is related to the arrow of time and the 2 nd law of of ther- 
modynamics p5| , occurs also in classical physics. How- 
ever, quantum mechanics produces a qualitatively new 
effect as well, known as decoherence jll],[36 37 1, sup- 
pressing off-diagonal elements in the reduced density ma- 
trices pi . This effect destroys the ability to observe long- 
range quantum superpositions within the subsystems, 
and is now rather well-understood and uncontroversial 
p0| , p8| -|42]| - the interested reader is referred to |l3| and 
a recent book on decoherence 0] for details. For in- 
stance, if our colloid was initially in a superposition of 
two locations separated by a centimeter, this macrosu- 
perposition would for all practical purposes be destroyed 
by the first collision with a water molecule, i.e., on a 
timescale Td oc of order r co ii, with the quantum superpo- 
sition surviving only on scales below the de de Broigle 
wavelength of the water molecules pHfhlJ] This means 



2 Decoherence picks out a preferred basis in the quantum- 
mechanical Hilbert space, termed the "pointer basis" by 
Zurek , in which superpositions are rapidly destroyed and 
classical behavior is approached. This normally includes the 
position basis, which is why we never experience superposi- 
tions of objects in macroscopically different positions. Deco- 
herence is quite generic. Although it has been claimed that 
this preferred basis consists of the maximal set of commuting 
observables that also commute with Hi nt (the "microstable 
basis" of Omnes ^3|), this is in fact merely a sufficient condi- 
tion, not a necessary one. If [Hint, a:] = for some observable 
x but [Hint,p] 7^ for its conjugate p, then the interaction 
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that Tdiss/idcc ~ M/m in our example, i.e., that decoher- 
ence is much faster than dissipation for macroscopic ob- 
jects, and this qualitative result has been shown to hold 
quite generally as well (see and references therein). 
Loosely speaking, this is because each microscopic par- 
ticle that scatters off of the subsystem carries away only 
a tiny fraction m/M of the total momentum, but essen- 
tially all of the necessary information. 
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FIG. 1. The qualitative behavior of a subsystem depends on 
the timescales for dynamics, dissipation and decoherence. This 
classification is by necessity quite crude, so the boundaries should 
not be thought of as sharp. 



C. Classification of systems 



by the effects of the surroundings. In the intermediate 
case where Td cc <C Td yn < Tdi ss , we have a familiar classi- 
cal system. 

The relation between t& cc and Tdi ss depends only on 
the form of flint, whereas the question of whether Td yn 
falls between these values depends on the normalization 
of Hi n t in equation (|J). Since Td oc ~ Tdi ss for microscopic 
(atom-sized) systems and Td cc "C Tdi SS for macroscopic 
ones, Figure 1 shows that whereas macroscopic systems 
can behave quantum-mechanically, microscopic ones can 
never behave classically. 



D. Three systems: subject, object and environment 

Most discussions of quantum statistical mechanics split 
the Universe into two subsystems |47j]: the object under 
consideration and everything else (referred to as the en- 
vironment). Since our purpose is to model the observer, 
we need to include a third subsystem as well, the subject. 
As illustrated in Figure 2, we therefore decompose the 
total system into three subsystems: 

• The subject consists of the degrees of freedom as- 
sociated with the subjective perceptions of the ob- 
server. This does not include any other degrees of 
freedom associated with the brain or other parts of 
the body. 

• The object consists of the degrees of freedom that 
the observer is interested in studying, e.g., the 
pointer position on a measurement apparatus. 

• The environment consists of everything else, i.e., 
all the degrees of freedom that the observer is not 
paying attention to. By definition, these are the 
degrees of freedom that we always perform a partial 
trace over. 



Let us define the dynamical timescale Td yn of a subsys- 
tem as that which is characteristic of its internal dynam- 
ics. For a planetary system or an atom, Td yn would be 
the orbital frequency. 

The qualitative behavior of a system depends on the 
ratio of these timescales, as illustrated in Figure 1. If 
Tdyn <C Tdcc, we are are dealing with a true quantum sys- 
tem, since its superpositions can persist long enough to 
be dynamically important. If Td yn 3> Tdi SS , it is hardly 
meaningful to view it as an independent system at all, 
since its internal forces are so week that they arc dwarfed 



will indeed cause decoherence for x as advertised. But this 
will happen even if [Hi n t,x] 7^ — all that matters is that 
[.ffiiit,j>] 7^ 0, i.e., that the interaction Hamiltonian contains 
("measures") x. 
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tion. It is the focus of the present paper, and the next 
section is devoted to quantitative calculations of decoher- 
ence in brain processes, aimed at determining whether 
the subject system should be classified as classical or 
quantum in the sense of Figure 1. We will return to 
Figure 2 and a more detailed discussion of its various 
subsystem interactions in Section [V. 



ENVIRONMENT 

H e 

(Always traced over) 



FIG. 2. An observer can always decompose the world into three 
subsystems: the degrees of freedom corresponding to her subjective 
perceptions (the subject), the degrees of freedom being studied (the 
object), and everything else (the environment). As indicated, the 
subsystem Hamiltonians H s , H a , H e and the interaction Hamilto- 
nians H BO , H oe , H Ee can cause qualitatively very different effects, 
which is why it is often useful to study them separately. This paper 
focuses on the interaction H se . 

Note that the first two definitions are very restrictive. 
Whereas the subject would include the entire body of 
the observer in the common way of speaking, only very 
few degrees of freedom qualify as our subject or object. 
For instance, if a physicist is observing a Stern-Gerlach 
apparatus, the vast majority of the ~ 10 28 degrees of 
freedom in the the observer and apparatus are counted 
as environment, not as subject or object. 

The term "perception" is used in a broad sense in item 
1, including thoughts, emotions and any other attributes 
of the subjectively perceived state of the observer. 

The practical usefulness in this decomposition lies in 
that one can often neglect everything except the object 
and its internal dynamics (given by H a ) to first order, 
using simple prescriptions to correct for the interactions 
with the subject and the environment. The effects of 
both H so and H oe have been extensively studied in the 
literature. H so involves quantum measurement, and gives 
rise to the usual interpretation of the diagonal elements of 
the object density matrix as probabilities. H oc produces 
decoherence, selecting a preferred basis and making the 
object act classically if the conditions in Figure I are met. 

In contrast, H sc , which causes decoherence directly in 
the subject system, has received relatively little atten- 



III. DECOHERENCE RATES 

In this section, we will make quantitative estimates 
of decoherence rates for neurological processes. We first 
analyze the process of neuron firing, widely assumed to be 
central to cognitive processes. We also analyze electrical 
excitations in microtubules, which Penrose and others 
have suggested may be relevant to conscious thought. 



A. Neuron firing 

Neurons (see Figure 3) are one of the key building 
blocks of the brain's information processing system. It is 
widely believed that the complex network of ~ 10 11 neu- 
rons with their nonlinear synaptic couplings is in some 
way linked to our subjective perceptions, i.e., to the sub- 
ject degrees of freedom. If this picture is correct, then if 
H s or H so puts the subject into a superposition of two 
distinct mental states, some neurons will be in a super- 
position of firing and not firing. How fast does such a 
superposition of a firing and non-firing neuron decohere? 

Let us consider this process in more detail. For in- 
troductory reviews of neuron dynamics, the reader is re- 
ferred to, e.g., [^8|-^0|. Like virtually all animal cells, 
neurons have ATP driven pumps in their membranes 
which push sodium ions out of the cell into the surround- 
ing fluids and potassium ions the other way. The former 
process is slightly more efficient, so the neuron contains a 
slight excess of negative charge in its "resting" state, cor- 
responding to a potential difference Uq ~ —0.07 V across 
the axon membrane ( "axolemma" ) . There is an inher- 
ent instability in the system, however. If the potential 
becomes substantially less negative, then voltage-gated 
sodium channels in the axon membrane open up, allow- 
ing Na + ions to come gushing in. This makes the poten- 
tial still less negative, causes still more opening, etc. This 
chain reaction, "firing" , propagates down the axon at a 
speed of up to 100 m/s, changing the potential difference 
to a value U\ that is typically of order +0.03 V ]49| . 

The axon quickly recovers. After less than ~ I ms, the 
sodium channels close regardless of the voltage, and large 
potassium channels (also voltage gated, but with a time 
delay) open up allowing K+ ions to flow out and restore 
the resting potential Uq. The ATP driven pumps quickly 
restore the Na + and K + concentrations to their initial 
values, making the neuron ready to fire again if triggered. 
Fast neurons can fire over I0 3 times per second. 
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FIG. 3. Schematic illustration of a neuron (left), a section of 
the myelinated axon (center) and and a piece of its axon membrane 
(right). The axon is typically insulated (myelinated) with small 
bare patches every 0.5 mm or so (so-called Nodes of Ranvier) where 
the voltage-sensitive sodium and potassium gates are concentrated 
|5l|j52| . If the neuron is in a superposition of firing and not firing, 
then N ~ 10 6 Na + ions are in a superposition of being inside and 
outside the cell (right). 

Consider a small patch of the membrane, assumed to 
be roughly flat with uniform thickness h as in Figure 3. 
If there is an excess surface density ±c of charge near 
the inside/outside membrane surfaces, giving a voltage 
differential U across the membrane, then application of 
Gauss' law tells us that a = eoE, where the electric field 
strength in the membrane is E = U/h and eo is the vac- 
uum permittivity. Consider an axon of length L and 
diameter d, with a fraction / of its surface area bare (not 
insulated with myelin). The total active surface area is 
thus A = ndLf, so the total number of Na + ions that 
migrate in during firing is 



N 



An = irdLfeojU! - U ) 
q qh 



(4) 



where q is the ionic charge (q — q e , the absolute value 
of the electron charge). Taking values typical for central 
nervous system axons p^ , p3|| , h — 8nm, d — 10 /im, 
L = 10cm, / = 10~ 3 , U = -0.07 V and Ux = +0.03 V 
gives N w 10 6 ions, and reasonable variations in our 
parameters can change this number by a few orders of 
magnitude. 



B. Neuron decoherence mechanisms 

Above we saw that a quantum superposition of the 
neuron states "resting" and "firing" involves of order a 
million ions being in a spatial superposition of inside and 
outside the axon membrane, separated by a distance of 
order h ~ 10 nm. In this subsection, we will compute the 
timescale on which decoherence destroys such a superpo- 
sition. 



In this analysis, the object is the neuron, and the su- 
perposition will be destroyed by any interaction with 
other (environment) degrees of freedom that is sensitive 
to where the ions are located. We will consider the fol- 
lowing three sources of decoherence for the ions: 

1. Collisions with other ions 

2. Collisions with water molecules 

3. Coloumb interactions with more distant ions 

There are many more decoherence mechanisms 10-^6) . 
Exotic candidates such as quantum gravity |54|| and 
modified quantum mechanics 55 5^] are generally much 
weaker [^6| . A number of decoherence effects may be even 
stronger than those listed, e.g., interactions as the ions 
penetrate the membrane — the listed effects will turn out 
to be so strong that we can make our argument by sim- 
ply using them as lower limits on the actual decoherence 
rate. 

Let p denote the density matrix for the position r of a 
single Na + ion. As reviewed in the Appendix, all three 
of the listed processes cause p to evolve as 



p(x, x', t + 1) = p(x, x', t )f(x, x', t) 



(5) 



for some function / that is independent of the ion state 
p and depends only on the interaction Hamiltonian H mt . 
This assumes that we can neglect the motion of the ion 
itself on the decoherence timescale — we will see that 
this condition is met with a broad margin. 



1. Ion-ion collisions 

For scattering of environment particles (processes 1 
and 2) that have a typical de Broigle wavelength A, we 
have |6| 



/(x,x',i) = e v 

( e -|x'-x| 2 Ai/2A 2 
~ 1 „-Ai 



for 
for 



x| < A, 
x| > A. 



(G) 



Here A is the scattering rate, given by A = n(av), where 
n is the density of scatterers, a is the scattering cross 
section and v is the velocity. The product av is aver- 
aged over a the velocity distribution, which we take to 
be a thermal (Boltzmann) distribution for correspond- 
ing to T = 37°Cw 310 K. The gist of equation (|) is 
that a single collision decoheres the ion down to the 
de Broigle wavelength of the scattering particle. The 
information I\2 communicated during the scattering is 
7i2 ~ log 2 (Ax/A) bits, where Ace is the initial spread in 
the position of our particle. 

Since the typical de Broigle wavelength of a Na + ion 
(mass m w 23m p ) or H 2 molecule (m « 18m p ) is 



5 



2irh 
\Z3mkT 



0.03 nm 



(7) 



at 310K, way smaller than the the membrane thickness 
h ~ 10 nm over which we need to maintain quantum 
coherence, we are clearly in the |x' — x| ^> A limit of 
equation (Q). This means that the spatial superposition 
of an ion decays exponentially A" 1 , of order its mean 
free time between collisions. Since the superposition of 
the neuron states "resting" and "firing" involves N such 
superposed ions, it thus gets destroyed on a timescale 

TEE (A^A)- 1 . 

Let us now evaluate r. Coulomb scattering between 
two ions of unit charge gives substantial deflection angles 
(0 ~ 1) with a cross section or orderf] 



(9) 



dipole as two opposing unit charges separated by a dis- 
tance y = p/q e <C b, so summing the two corresponding 
contributions from equation (0) gives a deflection angle 



2gq e P 
mv 2 b 2 ' 



This gives a cross section 



a = nb 



gqeP 



(ii) 



(12) 



for scattering with large (8 ~ 1) deflections. Although a 
is smaller than for the case of ion-ion collisions, n is larger 
because the concentration factor 77 drops out, giving a 
final result 



1 



V mkT 
Nnav Ngq e pn 



10 



-20 , 



(13) 



where v is the relative velocity and g ee l/47reo is the 
Coulomb constant. In thermal equilibrium, the kinetic 
energy mv 2 /2 is of order kT, so v ~ ^JkT/m. For the 
ion density, let us write n = nn^o, where the density 
of water molecules iih 2 o is about (1 g/cm 3 )/(18m p ) ~ 
10 23 /cm 3 and n is the relative concentration of ions (pos- 
itive and negative combined) . Typical ion concentrations 
during the resting state are [Na + ] =9.2 (120) mmol/1 and 
[K+] =140 (2.5) mmol/1 inside (outside) the axon mem- 
brane |Q , corresponding to total Na + + K + concentra- 
tions of 77 w 0.00027 (0.00022) inside (outside). To be 
conservative, we will simply use r\ sa 0.0002 throughout. 
Ion-ion collisions therefore destroy the superposition on 
a timescale 



1 



Nnav 



\J m(kT) s 
Ng 2 q^n 



10- 20 s. 



(10) 



2. Ion-water collisions 

Since H2O molecules are electrically neutral, the cross- 
section is dominated by their electric dipole moment 
p w 1.85Debye« (0.0385 nm) x q e . We can model this 



3. Interactions with distant ions 

As shown in the Appendix, long-range interaction with 
a distant (environment) particle gives 



f{r,r',t)=p 2 [M(r'-r)t/h] 



(14) 



up to a phase factor that is irrelevant for decoherence. 
Here p 2 is the Fourier transform of pi (r) ee p 2 (r, r) , the 
probability distribution for the location of the environ- 
ment particle. M is the 3x3 Hessian matrix of second 
derivatives of the interaction potential of the two parti- 
cles at their mean separation. A slightly less general for- 
mula was derived in the seminal paper f45|| . For roughly 
thermal states, P2 (and thus p) is likely to be well ap- 
proximated by a Gaussian |5^p^| . This gives 



/(r,r',t) = e -5( r '- r ) tMtSM ( r '- r )* 2 /* 2 



(15) 



where X = (^2^2) ~ ( r ^)( r 2) ^ s the covariance matrix of 
the location of the environment particle. Decoherence 
is destroyed when the exponent becomes of order unity, 
i.e., on a timescale 



[(r' - r)'M'53M(r' - r)] 1/2 ft. 



(16) 



3 If the first ion starts at rest at ri = (0, 0, 0) and the sec- 
ond is incident with T2 = (vt, b, 0), then a very weak scatter- 
ing with deflection angle 9 < 1 will leave these trajectories 
roughly unchanged, the radial force F = gq 2 /|ri — r2| 2 merely 
causing a net transverse acceleration |57j 



AVy = 



-dt 



gq 2 bdt 



[b 2 + (vt) 2 ]*/ 2 mvb' 



(«) 



The approximation breaks down as the deflection angle 9 w 
Avy/v approaches unity. This occurs for b ~ gq 2 /mv 2 , giving 
a — Tvb 2 as in equation (^). 



Assuming a Coulomb potential V = gq 2 j\v2 — Y\\ gives 



M 



(3aa* — V)gq- jar where a ee V2 — v\ 



1. For thermal states, we have the isotropic case S 
(Aie) 2 I, so equation (|TJ) reduces to 



ha 3 



gq 2 \r' 



r\Ax 



(1 + 3 cos 2 e) 



-1/2 



(17) 



where cos6> ee a • (r' — r)/|r' — r|. To be conservative, 
we take Ax to be as small as the uncertainty principle 
allows. With the thermal constraint (Ap) 2 /m < kT on 
the momentum uncertainty, this gives 



G 



Ax = 



2Ap yfmkT' 



(18) 



Substituting this into equation ( |l7j ) and dividing by the 
number of ions N , we obtain the dccohcrcncc timescale 



^\/mkT 



Ngq 2 



(19) 



caused by a single environment ion a distance a away. 
Each such ion will produce its own suppression factor /, 
so we need to sum the exponent in equation ([l5]) over all 
ions. Since the tidal force M oc a~ 3 causes the exponent 
to drop as a -6 , this sum will generally be dominated by 
the very closest ion, which will typically be a distance 
a ~ n -1 / 3 away. We are interested in decoherence for 
separations |r' — r| = h, the membrane thickness, which 
gives 



ymfcT 
Ngq 2 nh 



10" 19 s. 



(20) 



The relation between these different estimates is dis- 
cussed in more detail in the Appendix. 



C. Microtubules 

Microtubules are a major component of the cytoskele- 
ton, the "scaffolding" that helps cells maintain their 
shapes. They are hollow cylinders of diameter D = 
24 nm made up of 13 filaments that are strung together 
out of proteins known as tubulin dimers. These dimers 
can make transitions between two states known as a 
and 0, corresponding to different electric dipole moments 
along the axis of the tube. It has been argued that micro- 
tubules may have additional functions as well, serving as 
a means of energy and information transfer pC| ] . A model 
has been presented whereby the dipole-dipole interac- 
tions between nearby dimers can lead to long-range po- 
larization and kink-like excitations that may travel down 
the microtubules at speeds exceeding 1 m/s [ |iio| . 

Penrose has gone further and suggested that the dy- 
namics of such excitations can make a microtubule act 
like a quantum computer, and that microtubules are the 
site of of human consciousness (|] . This idea has been fur- 
ther elaborated |H] [24|] employing methods from string 
theory, with the conclusion that quantum superpositions 
of coherent excitations can persist for as long as a second 
before being destroyed by decoherence. See also |6l|,|62[ . 
This was hailed as a success for the model, the interpre- 
tation being that the quantum gravity effect on micro- 
tubules was identified with the human though process on 
this same timescale. 

This decoherence rate r ~ 1 s was computed assuming 
that quantum gravity is the main decoherence source. 
Since this quantum gravity model is somewhat contro- 
versial |32| and its effect has been found to be more than 



20 orders of magnitude weaker than other decoherence 
sources in some cases ft6f , it seems prudent to evalu- 
ate other decoherence sources for the microtubule case 
as well, to see whether they are in fact dominant. We 
will now do so. 

Using coordinates where the x-eods is along the tube 
axis, the above-mentioned models all focus on the time- 
evolution of p(x), the average x-component of the electric 
dipole moment of the tubulin dimers at each x. In terms 
of this polarization function p(x), the net charge per unit 
length of tube is — p'(x). The propagating kink- like exci- 
tations 16011 are of the form 



p(x) 



+po for x <C x , 
—p for x > xq, 



(21) 



where the kink location xq propagates with constant 
speed and has a width of order a few tubulin dimers. 
The polarization strength po is such that the total charge 
around the kink is Q = — J p'{x)dx = 2po ~ 940g e , due 
to the presence of 18 Ca 2+ ions on each of the 13 fila- 
ments contributing to po |30| . 

Suppose that such a kink is in two different places 
in superposition, separated by some distance |r' — r|. 
How rapidly will the superposition be destroyed by de- 
coherence? To be conservative, we will ignore colli- 
sions between polarized tubulin dimers and nearby water 
molecules, since it has been argued that these may be in 
some sense ordered and part of the quantum system ]24| ] 
- although this argument is difficult to maintain for the 
water outside the microtubule, which permeates the en- 
tire cell volume. Let us instead apply equation (|l^), with 
N = Q/q e ~ 10 3 . The distance to the nearest ion will 
generally be less than a = R + n^ 1 / 3 ~ 26 nm, where the 
tubulin diameter D = 24 nm dominates over the inter- 
ion separation n -1 / 3 ~ 2nm in the fluid surrounding 
the microtubule. Superpositions spanning many tubuline 
dimers (|r' — r| ^> D) therefore decohere on a timescale 
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-13 , 



(22) 



due to the nearest ion alone. This is quite a conserva- 
tive estimate, since the other nD 3 ~ 10 3 ions that are 
merely a small fraction further away will also contribute 
to the decoherence rate, but it is nonetheless 6-7 orders 
of magnitude shorter than the estimates of Mavromatos 
& Nanopoulos p5|-|27||. We will comment on screening 
effects below. 



1. Decoherence summary 

Our decoherence rates are summarized in Table 1 . How 
accurate are they likely to be? 

In the calculations above, we generally tried to be con- 
servative, erring on the side of underestimating the deco- 
herence rate. For instance, we neglected that N potas- 
sium ions also end up in superposition once the neuron 
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firing is quenched, we neglected the contribution of other 
abundant ions such as CI - to rj, and and we ignored col- 
lisions with water molecules in the microtubule case. 

Since we were only interested in order-of-magnitudc 
estimates, we made a number of crude approximations, 
e.g., for the cross sections. We neglected screening ef- 
fects because the decoherence rates were dominated by 
the particles closest to the system, i.e., the very same par- 
ticles that are responsible for screening the charge from 
more distant ones. 



Table 1. Decoherence timescales. 



Object 


Environment 


^dcc 


Neuron 


Colliding ion 


10" 20 s 


Neuron 


Colliding H 2 


10- 20 s 


Neuron 


Nearby ion 


10" 19 s 


Microtubule 


Distant ion 


icr 13 s 



to the limits of plausibility, it is clearly impossible to 
change the basic conclusion that Td cc "C 10~ 3 s, i.e., that 
we are dealing with a classical system in the sense of Fig- 
ure 1. Even the tiniest neuron imaginable, with only a 
single ion (N = 1) traversing the cell wall during firing, 
would have Td cc ~ 10~ 14 s. Likewise, reducing the ef- 
fective microtubule kink charge to a small fraction of q e 
would not help. 

How are we to understand the above-mentioned claims 
that brain subsystems can be sufficiently isolated to 
exhibit macroquantum behavior? It appears that the 
subtle distinction between dissipation and decoherence 
timescales has not always been appreciated. 

B. Implications for the subject-object-environment 
decomposition 



IV. DISCUSSION 

A. The classical nature of brain processes 

The calculations above enable us to address the ques- 
tion of whether cognitive processes in the brain consti- 
tute a classical or quantum system in the sense of Fig- 
ure 1. If we take the characteristic dynamical timescale 
for such processes to be Td yn ~ 10~ 2 s — 10° s (the ap- 
parent timescale of e.g., speech, thought and motor re- 
sponse) , then a comparison of Td yn with Td cc from Table 1 
shows that processes associated with either conventional 
neuron firing or with polarization excitations in micro- 
tubules fall squarely in the classical category, by a mar- 
gin exceeding ten orders of magnitude. Neuron firing it- 
self is also highly classical, since it occurs on a timescale 
Tdy n ~ 10~ 3 — 10~ 4 s |53|| . Even a kink- like microtubule 
excitation is classical by many orders of magnitude, since 
it traverses a short tubule on a timescale Td yn ~ 5 x 10~ 7 s 
& 

What about other mechanisms? It is worth noting 
that if (as is commonly believed) different neuron fir- 
ing patterns correspond in some way to different con- 
scious perceptions, then consciousness itself cannot be 
of a quantum nature even if there is a yet undiscovered 
physical process in the brain with a very long decoherence 
time. As mentioned above, suggestions for such candi- 
dates have involved, e.g., superconductivity Jl2| , super- 
fluidity [|l3| , electromagnetic fields flif| , Bose condensa- 
tion | 1 5 U l6|l, superflourescence fl7| and other mechanisms 
^8| , |l9| . The reason is that as soon as such a quantum 
subsystem communicates with the constantly decohering 
neurons to create conscious experience, everything deco- 
heres. 

How extreme variations in the decoherence rates can 
we obtain by changing our model assumptions? Although 
the rates can be altered by a few of orders of magnitudes 
by pushing parameters such as the neuron dimensions, 
the myelination fraction or the microtubule kink charge 



Let us now discuss the subsystem decomposition of 
Figure 2 in more detail in light of our results. As the 
figure indicates, the virtue of this decomposition into 
subject, object and environment is that the subsystem 
Hamiltonians H s , H a , H e and the interaction Hamiltoni- 
ans H so , H oc , H se can cause qualitatively very different 
effects. Let us now briefly discuss each of them in turn. 

Most of these processes are schematically illustrated 
in Figure 4 and Figure 5, where for purposes of illus- 
tration, we have shown the extremely simple case where 
both the subject and object have only a single degree of 
freedom that can take on only a few distinct values (3 
for the subject, 2 for the object). For definiteness, we 
denote the three subject states | "), and |^), and in- 
terpret them as the observer feeling neutral, happy and 
sad, respectively. We denote the two object states | f ) 
and ||), and interpret them as the spin component ("up" 
or "down") in the z-direction of a spin- 1/2 system, say a 
silver atom. The joint system consisting of subject and 
object therefore has only 2x3 = 6 basis states: |" |), 
| I'D, |ot), |^T), |^4). In Figures 4 and 5, we 

have therefore plotted p as a 6 x 6 matrix consisting of 
nine two- by- two blocks. 
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FIG. 4. Time evolution of the 6x6 density matrix for the basis 
states | !"["), of), |^T)i 1^1) as the object evolves in 

isolation, then decoheres, then gets observed by the subject. The 
final result is a statistical mixture of the states |oT) an d 
simple zero-entropy states like the one we started with. 



1 . Effect of H : constant entropy 

If the object were to evolve during a time interval t 
without interacting with the subject or the environment 
(H so = H oc = 0), then according to equation (nj) its 
reduced density matrix p Q would evolve into U p U with 
the same entropy, since the time-evolution operator U = 
e -iH„t j g unwary. 

Suppose the subject stays in the state 1 1' ) and the 
object starts out in the pure state ||). Let the object 
Hamiltonian H Q correspond to a magnetic field in the y- 
direction causing the spin to precess to the x-direction, 
i.e., to the state (|T) + li))/v / 2- The object density matrix 
p then evolves into 

1 



Po = u\m^ = 2(lt) + l4))«TI + ai) 

= i(iT>ai + it>ai + ii>ai + ii>ai) ) 



(23) 



corresponding to the four entries of 1/2 in the second 
matrix of Figure 4. 

This is quite typical of pure quantum time evolution: a 
basis state eventually evolves into a superposition of ba- 
sis states, and the quantum nature of this superposition 
is manifested by off-diagonal elements in p Q . Another fa- 
miliar example of this is the familiar spreading out of the 
wave packet of a free particle. 



Figure 4 illustrates the case of an ideal measurement, 
where the subject starts out in the state | Z) and H so is of 
such a form that gets perfectly correlated with the object. 
In the language of Section ||, an ideal measurement is a 
type of communication where the mutual information I12 
between the subject and object systems is increased to its 
maximum possible value. Suppose that the measurement 
is caused by _ff so becoming large during a time interval so 
brief that we can neglect the effects of H s and H Q . The 
joint subject+object density matrix p so then evolves as 
p so t-y Up so U', where U = exp [— i J H so dt] . If observing 
|t) makes the subject happy and ||) makes the subject 
sad, then we have U\ 1*T) =JoT) and U\ I'D = The 
state given by equation would therefore evolve into 



A, = ^(|:>(:|)®(|T><TI + li>ai)^ 



= -(u\:V(:^ + u\:i)(:i\u^ 



= -(|ot)(ot| + \M)(M 



(25) 
(26) 
(27) 



as illustrated in Figure 4. This final state contains a 
mixture of two subjects, corresponding to definite but 
opposite knowledge of the object state. According to 
both of them, the entropy of the object has decreased 
from one bit to zero bits. 

In general, we see that the object decreases its en- 
tropy when it exchanges information with the subject 
and increases when it exchanges information with the 
environment.^ Loosely speaking, the entropy of an ob- 
ject decreases while you look at it and increases while 
you don'tR. 



2. Effect of H oa : increasing entropy 

This was the effect of H Q alone. In contrast, H oe will 
generally cause decoherence and increase the entropy of 
the object. As discussed in detail in Section III and the 
Appendix, it entangles it with the environment, which 
suppresses the off-diagonal elements of the reduced den- 
sity matrix of the object as illustrated in Figure 4. If H oc 
couples to the z-component of the spin, this destroys the 
terms |t)(l| and ||)(t|- Complete decoherence therefore 
converts the final state of equation (E31) into 



Po = i(ITXTi + UX1I), 



(24) 



corresponding to the two entries of 1/2 in the third ma- 
trix of Figure 4. 



3. Effect of H so : decreasing entropy 

Whereas H oe typically causes the apparent entropy of 
the object to increase, H so typically causes it to decrease. 



4 If n bits of information are exchanged with the environ- 
ment, then equation (^) shows that the object entropy will 
increase by this same amount if the environment is in ther- 
mal equilibrium (with maximal entropy) throughout. If we 
were to know the state of the environment initially (by our 
definition of environment, we do not), then both the object 
and environment entropy will typically increase by n/2 bits. 

5 Here and throughout, we are assuming that the total 
system, which is by definition isolated, evolves according to 
the Schrodinger equation (|l|). Although modifications of the 
Schrodinger equation have been suggested by some authors, 
either in a mathematically explicit form as in J5f|^(| or ver- 
bally as a so-called reduction postulate, there is so far no 
experimental evidence suggesting that modifications are nec- 
essary. The original motivations for such modifications were 

1. to be able to interpret the diagonal elements of the 
density matrix as probabilities and 

2. to suppress off-diagonal elements of the density matrix. 

The subsequent discovery by Everett ^iij that the probability 
interpretation automatically appears to hold for almost all 
observers in the final superposition solved problem 1, and is 
discussed in more detail in, e.g., pj|,|66|-[74 1 . The still more 
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FIG. 5. Time evolution of the same 6x6 density matrix as in 
Figure 4 when the subject evolves in isolation, then decoheres. The 
object remains in the state |j") the whole time. The final result is 
a statistical mixture of the two states |oT) an d I^T)- 



4- Effect if H s : the thought process 

So far, we have focused on the object and discussed 
effects of its internal dynamics (H a ) and its interactions 
with the environment (H oc ) and subject (H so ). Let us 
now turn to the subject and consider the role played by 
its internal dynamics (H s ) and interactions with the en- 
vironment (H se ). In his seminal 1993 book, Stapp H 
presents an argument about brain dynamics that can be 
summarized as follows. 

1. Since the brain contains ~ 10 11 synapses connected 
together by neurons in a highly nonlinear fashion, 
there must be a huge number of metastable rever- 
berating patters of pulses into which the brain can 
evolve. 

2. Neural network simulations have indicated that the 
metastable state into which a brain does in fact 
evolves depends sensitively on the initial conditions 
in small numbers of synapses. 

3. The latter depends on the locations of a small num- 
ber of calcium atoms, which might be expected to 
be in quantum superpositions. 

4. Therefore, one would expect the brain to evolve 
into a quantum superposition of many such 
metastable configurations. 

5. Moreover, the fatigue characteristics of the synap- 
tic junctions will cause any given metastable state 



to become, after a short time, unstable: the 
subject will then be forced to search for a new 
metastable configuration, and will therefore con- 
tinue to evolve into a superposition of increasingly 
disparate states. 

If different states (perceptions) of the subject correspond 
to different metastable states of neuron firing patterns, a 
definite perception would eventually evolve into a super- 
position of several subjectively distinguishable percep- 
tions. 

We will follow Stapp in making this assumption about 
H s . For illustrative purposes, let us assume that this can 
happen even at the level of a single thought or snap de- 
cision where the outcome feels unpredictable to us. Con- 
sider the following experiment: the subject starts out 
with a blank face and counts silently to three, then makes 
a snap decision on whether to smile or frown. The time- 
evolution operator U = exp \—i J H s dt\ will then have 

the property that U\Z) — (|c) + \^))/y/2, so the sub- 
ject density matrix p s will evolve into 

Ps = u\:){:\tf = ~(\c) + \^))({^\ + {M) 



= -(|o)<o| + |, 



+ + (28) 



corresponding to the four entries of 1/2 in the second 
matrix in Figure 5. 



5. Effect of H Be : subject decoherence 

Just as H oc can decohere the object, H se can decohere 
the subject. The difference is that whereas the object can 
be either a quantum system (with small H oc ) or a classi- 
cal system (with large H oe ), a human subject always has 
a large interaction with the environment. As we showed 
in Section III, Td oc <C Td yn for the subject, 
feet of H, 



Tdoc "C Tdyn tor trie subject, i.e., the ef- 
is faster than that of H s by many orders of 
magnitude. This means that we should strictly spea king 
not think of macrosuperpositions such as equation (|2q) 
as first forming and then decohering as in Figure 5 - 
rather, subject decoherence is so fast that such superpo- 
sitions decohere already during their process of forma- 
tion. Therefore we are never even close to being able to 
perceive superpositions of different perceptions. Reduc- 
ing object decoherence (from H oc ) during measurement 
would make no difference, since decoherence would take 
place in the brain long before the transmission of the ap- 
propriate sensory input through sensory nerves had been 
completed. 



recent discovery of decoherence ]1 l| , |36|j37[ ] solved problem 2, C. H e and H aoe 
as well as explaining so-called superselection rules for the first 

time (why for instance the position basis has a special status) mi • j. ■ r ,, , , , 

.r-| v ' I ne environment is ot course the most complicated sys- 

'' tern, since it contains the vast majority of the degrees of 
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freedom in the total system. It is therefore very fortu- 
nate that we can so often ignore it, considering only those 
limited aspects of it that affect the subject and object. 

For the most general H, there can also be an ugly 
irreducible residual term H soe = H — H s — H — H c — 

H-so ^oc H se . 



D. Implications for modeling cognitive processes 

For the neural network community, the implication of 
our result is "business as usual", i.e., there is no need 
to worry about the fact that current simulations do not 
incorporate effects of quantum coherence. The only rem- 
nant from quantum mechanics is the apparent random- 
ness that we subjectively perceive every time the subject 
system evolves into a superposition as in equation (28), 
but this can be simply modeled by including a random 
number generator in the simulation. In other words, the 
recipe used to prescribe when a given neuron should fire 
and how synaptic coupling strengths should be updated 
may have to involve some classical randomness to cor- 
rectly mimic the behavior of the brain. 



1. Hyper- classicality 

If a subject system is to be a good model of us, H so 
and H so need to meet certain criteria: decoherence and 
communication are necessary, but fluctuation and dissi- 
pation must be kept low enough that the subject does 
not lose its autonomy completely. 

In our study of neural processes, we concluded that the 
subject is not a quantum system, since Td oc <C Td yn - How- 
ever, since the dissipation time Tdi ss for neuron firing is 
of the same order as its dynamical timescale, we see that 
in the sense of Figure 1, the subject is not a simple clas- 
sical system either. It is therefore somewhat misleading 
to think of it as simply some classical degrees of freedom 
evolving fairly undisturbed (only interacting enough to 
stay decohered and occasionally communicate with the 
outside world). Rather, the semi-autonomous degrees of 
freedom that constitute the subject are to be found at a 
higher level of complexity, perhaps as metastable global 
patters of neuron firing. 

These degrees of freedom might be termed "hyper- 
classical" : although there is nothing quantum- 
mechanical about their equations of motion (except that 
they can be stochastic) , they may bear little resemblance 
with the underlying classical equations from which they 
were derived. Energy conservation and other familiar 
concepts from Hamiltonian dynamics will be irrelevant 
for these more abstract equations, since neurons are en- 
ergy pumped and highly dissipative. Other examples of 
such hyper-classical systems include the time-evolution 
of the memory contents of a regular (highly dissipative) 



digital computer as well as the motion on the screen of 
objects in a computer game. 



2. Nature of the subject system 

In this paper, we have tacitly assumed that conscious- 
ness is synonymous with certain brain processes. This is 
what Lockwood terms the "identity theory" |3(| . It dates 
back to Hobbes (^1660) and has been espoused by, e.g., 
Russell, Feigl, Smart, Armstrong, Churchland and Lock- 
wood himself. Let us briefly explore the more specific 
assumption that the subject degrees of freedom are our 
perceptions. In this picture, some of the subject degrees 
of freedom would have to constitute a "world model", 
with the interaction H so such that the resulting commu- 
nication keeps these degrees of freedom highly correlated 
with selected properties of the outside world (object + 
environment). Some such properties, i.e., 

• the intensity of the electromagnetic on the retina, 
averaged through three narrow-band filters (color 
vision) and one broad-band filter (black-and-white 
vision) , 

• the spectrum of air pressure fluctuations in the ears 
(sound), 

• the chemical composition of gas in the nose (smell) 
and solutions in the mouth (taste), 

• heat and pressure at a variety of skin locations, 

• locations of body parts, 

are tracked rather continuously, with the corresponding 
mutual information 7 12 between subject and surround- 
ings remaining fairly constant. Persisting correlations 
with properties of the past state of the surroundings 
(memories) further contribute to the mutual information 
Ii2- Much of I12 is due to correlations with quite subtle 
aspects of the surroundings, e.g., the contents of books. 
The total mutual information I\2 between a person and 
the external world is fairly low at birth, gradually grows 
through learning, and falls when we forget. In contrast, 
most innate objects have a very small mutual informa- 
tion with the rest of the world, books and diskettes being 
notable exceptions. 

The extremely limited selection of properties that the 
subject correlates with has presumably been determined 
by evolutionary utility, since it is known to differ between 
species: birds perceive four primary colors but cats only 
one, bees perceive light polarization, etc. In this picture, 
we should therefore not consider these particular ( "classi- 
cal" ) aspects of our surroundings to be more fundamental 
than the vast majority that the subject system is uncor- 
rected with. Morover, our perception of e.g. space is as 
subjective as our perception of color, just as suggested 
by e.g. M. 
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3. The binding problem 

One of the motivations for models with quantum co- 
herence in the brain was the so-called binding problem. 
In the words of James |75|,|76| , "the only realities are the 
separate molecules, or at most cells. Their aggregation 
into a 'brain' is a fiction of popular speech" . James' con- 
cern, shared by many after him, was that consciousness 
did not seem to be spatially localized to any one small 
part of the brain, yet subjectively feels like a coherent 
entity. Because of this, Stapp || and many others have 
appealed to quantum coherence, arguing that this could 
make consciousness a holistic effect involving the brain 
as a whole. 

However, non-local degrees of freedom can be impor- 
tant even in classical physics, For instance, oscillations 
in a guitar string are local in Fourier space, not in real 
space, so in this case the "binding problem" can be solved 
by a simple change of variables. As Eddington remarked 
|77|l , when observing the ocean we perceive the moving 
waves as objects in their own right because they display a 
certain permanence, even though the water itself is only 
bobbing up and down. Similarly, thoughts are presum- 
ably highly non-local excitation patterns in the neural 
network of our brain, except of a non-linear and much 
more complex nature. In short, this author feels that 
there is no binding problem. 

4- Outlook 

In summary, our decoherence calculations have in- 
dicated that there is nothing fundamentally quantum- 
mechanical about cognitive processes in the brain, sup- 
porting the Hepp's conjecture |3^| . Specifically, the com- 
putations in the brain appear to be of a classical rather 
than quantum nature, and the argument by Lisewski f7§| l 
that quantum corrections may be needed for accurate 
modeling of some details, e.g., non-Markovian noise in 
neurons, does of course not change this conclusion. This 
means that although the current state-of-the-art in neu- 
ral network hardware is clearly still very far from be- 
ing able to model and understand cognitive processes as 
complex as those in the brain, there are no quantum me- 
chanical reasons to doubt that this research is on the 
right track. 
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APPENDIX: DECOHERENCE FORMULAS 

The quantitative effect of decoherence from both short 
range interactions (scattering) and long-range interac- 
tions was first derived in a seminal paper by Joos & Zeh 
p5[ . Since our application involved scattering between 
particles of comparable mass, we used a generalized ver- 
sion of these results that included the effect of recoil jl6) . 
In this Appendix, we derive a slightly generalized formula 
for long-range interactions, and briefly comment on the 
relation between these short-range and long-range limit- 
ing cases. 

1. Decoherence due to tidal forces 

Even if the dissipation and fluctuation caused by H lnt 
is dynamically unimportant, Hi and H 2 can be neglected 
in equation when calculating the decoherence effect in 
the many cases where the interaction Hamiltonian deco- 
heres the object on a timescale far below the dynamical 
time. In this approximation, we consider two particles 
with an interaction H = H mt = V(j2 — for some 
potential V. According to equation p), the two-particle 
density matrix p therefore evolves as 

p(r 1 ,r' 1 ,r 2 ,r' 2 ,tQ + t) 
= p(n, ri, r 2 , r' 2 , t) e - i [ v te- r O-i'(*' a -«i)]/\ (Al) 

Following [H , we assume that the two particles are fairly 
localized near their initial average positions 

r? = <r i ) =tr[r < p j (t )] > (A2) 

i = 1,2, and approximate the potential by its second 
order Taylor expansion 

U(r 2 - n) w V(a) - F • (x 2 - Xi) 

+ i(x 2 -x 1 ) t M(x 2 -x 1 ). (A3) 

Here F == — VU(a) is the average force, M is the Hes- 
sian matrix My = didjV(a) and a = r° — rj. We have in- 
troduced relative coordinates x; = r^— r?. Assuming that 
the two particles are independent initially as in [^9], i.e., 
that p(to) takes the separable form p(xi, x^, x 2 , x 2 , to) — 
pi(xi,x' 1 ,to)P2(x2,x 2 ,t ), this gives 

/^(x^x'^io +t) = tr 2( o(t +t) = 

/"pfxijx'^x,^ +t)d 3 x = / 9 1 (x 1 ,x' 1) i )/(xi,x' 1 ,t), (A4) 
where 
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/(xi,xi,i) « 

e Wu<,t) J (O2 (x 2 ,x' 2 ,t ) e - lt ( x 'i- Xl ) tMx2 /' i d 3 x 2 = 
e^ (xi ' x 'i'* ) p 2 [M(x' 1 - Xi)t/fc]. (A5) 
Here the phase factor 

e i0(x,x',t) E e f [F(x'-x)+Ix"Mx'-ix'Mx] ,£gs 

is of no importance for decohcrcncc, since it does not 
suppress the magnitude |pi(xi, x^, t)\ of the off-diagonal 
elements - it merely causes momentum transfer related 
to fluctuation and dissipation. It is the other term 
that causes decoherence. p 2 is the Fourier transform of 
p 2 (x) = p 2 (x, x, to), the probability distribution for the 
location of the environment particle. 

2. Properties of the effect 

Let us briefly discuss some qualitative features of equa- 
tion ( |A"5| ). Since p 2 (0) = J p 2 (x 2 )d 3 a; 2 = tr,o 2 = 1, 
Pi(x, x') remains unchanged on the diagonal x = x'. 
This is because H lrLt is not changing the position of our 
our object particle, merely its momentum. Since the 
mean position (x 2 ) = f p 2 x 2 <j 3 z 2 — tr [x 2 p 2 ] = van- 
ishes (using equation (A2)), we have Vp 2 (0) = 0. In 
fact, |/| takes a maximum on the diagonal, and the 
Riemann-Lebesgue Lemma shows that |/| = |p 2 | < 1 
whenever x ^ x', with equality only for the unphys- 
ical case where pi is a delta function, i.e., where the 
location of the environment particle is perfectly known. 
^6^1/(0)1 = -M<x 2 x 2 )Mi 2 /2?i 2 , so so the larger (x 2 x 2 ) 
is (i.e., the more spread out the environment particle is), 
the closer to the diagonal decoherence will suppress our 
density matrix. 

Since M is the shear matrix of the force field — W, we 
see that it is tidal forces that are causing the decoherence 
- the average force F simply contributes to the phase 
factor e 1 ^ . Specifically, the rate at which our object de- 
grees of freedom ri decohere grows with the tidal force 
that it exerts on the environment: if the environment 
particle is spread out with (x 2 x 2 ) large, experiencing a 
wide range of forces from the object, object decoherence 
is rapid. In the opposite situation, where the object is 
spread out and the environment is not, the object will 
experience strong classical tidal forces but no decoher- 
ence. 



3. Relation between long-range and short-range 
decoherence 



Above we derived the effect of decoherence from long- 
range tidal forces. Another interesting case that has been 
solved analytically Q is that of short-range interactions 



that can be modeled as scattering events. If the scatter- 
ing takes place during short enough a time interval that 
we can neglect the internal dynamics of the object, then 
its reduced density matrix changes as ||46| 



/9i(r,r') i ► pi(t,t')p 



r — r 



(A7) 



where p(q) is the probability distribution for the momen- 
tum transfer q in the collision. This equation generalizes 
the scattering result of |||] by including the effect of re- 
coil. The larger the uncertainty in momentum transfer, 
the stronger the decoherence effect becomes, since widen- 
ing p narrows its Fourier transform p. Changing the mean 
momentum transfer (q) does not affect the decoherence, 
merely contributes a phase factor j ust as F did above. 
Typically, the last factor in equation (A7) destroys coher- 
ence down to scales of order the de Broigle wavelength 
of the scatterer, with directional modulations from the 
angular dependence of the scattering cross section. Gen- 
eralization to a steady flux of scattering particles ]46| ] 
gives equatio n (p| ). 

Equation ( jA7| ) has str ikin g similarities with the tidal 
force result of equation (A5): in both cases, the density 
matrix gets multiplied by the Fourier transform of a prob- 
ability distribution. If fact, up to u nin teresting phase 
factors, we can re writ e our equation (A5) in exactly the 
form of equation ( A7) by redefining p to be the probabil- 
ity distribution for momentum transfer q = M(x 2 — xi)t 
due to tidal forces for a fixed Xi, i.e., 



P(q) =P2(x 2 ) 



d 3 x 2 _ p 2 (xi+M- 1 q/t) 
d?q ~ i 3 detM 



(A8) 



Fourier transforming this expression and substituting the 
result into equation (A7), we recover equation ( |A5| ) up 
to a phase factor. 

Perhaps the simplest way to understand all these re- 
sults is in terms of Wigner functions |7{| . If PF(x! , p x ) is 
the Wigner phase space distribution for the object parti- 
cle, then any of the momentum-transferring interactions 
that we have considered will take the form 



PF(xi,pi)^ J IF(xi,pi - q)p(q,xi)d 3 q 



(A9) 



for some probability distribution p that may or may not 
depend on Xi . Since the density matrix 

Pi(xi,x' 1 )= Jw(^±^, P y-^-^d 3 p (AlO) 

is just the Wigner function Fourier transformed in the 
momentum direction (an d ro tated by 45°), the convolu- 
tion with p in equation ( |A9| ) reduces to a simple multi- 
plication with p in equation (A7). 
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